Orbital symmetry fingerprints for magnetic adatoms in graphene 
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In this paper, we describe the formation of local resonances in graphene in the presence of magnetic 
adatoms containing localized orbitals of arbitrary symmetry, corresponding to any given angular 
momentum state. We show that quantum interference effects which are naturally inbuilt in the 
honeycomb lattice in combination with the specific orbital symmetry of the localized state lead 
to the formation of distinctive fingerprints in the local density of states and in the differential 
conductance curves, allowing scanning tunneling probes to characterize adatoms and defects in 
graphene. 
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INTRODUCTION 



Graphene is a single atomic layer of graphite whose 
low energy quasiparticles behave as massless Dirac 
fermions^^. As an open surface, graphene offers a solid 
playground for the detection and local manipulation of 
quantum states with scanning tunneling (STM) probes. 
This perspective is particularly promising for adatoms, 
which can be dragged with atomic precision^ and can 
have their magnetic state monitored and controlled with 
the application of an external gate voltage^^^. There has 
been substantial progress in the quality of the STM ex- 
periments in graphene in the last few years^"— . Recent 
experiments reported the observation of Landau levels 
spontaneously generated by strain on the top of nanobub- 
bles in grapheneii, and the observation of charge polar- 
ization effects around a Co adatom^^. 

Although the microscopic theory of STM is well under- 
stood in metallic hostsi^^^U, in graphene the sublattice 
quantum numbers play a role in the interference effects 
that drive the emergence of Fano resonances^ ^ nearby 
adatoms, in the presence of an STM tip. In particu- 
lar, for adatoms that sit at the center of the honeycomb 
hexagon {H site), destructive interference between the 
different electronic paths of hybridization with the two 
sublattices may give rise to a suppression of the Fano 
resonance of the localized state^^Ti^i, and also change the 
scattering rate of the localized electrons due to the pres- 
ence of the fermionic bath^^. In general, the broadening 
of a localized state, magnetic or not, is expected to scale 
as A{uj) oc |cj|^, where r is the scaling dimension of the 
DOS of the host material, which in graphene is r = 1 
(r = for metals). In graphene, nevertheless, localized 
orbitals located either in substitutional impurity sites {S 
sites) or in H sites and which also preserve the Csy point 
group symmetry of each sublattice are effectively damped 
at low energies by a fermionic bath with r = 322,25^ 



to quantum interference effects. This effect suggests that 
the local density of states (LDOS) can by quite suscepti- 
ble to the orbital symmetry of the localized state, allow- 
ing STM probes to characterize adatoms and defects in 
graphene. 

In this work, we describe in detail the effect of the 
localized orbital symmetry into the emergence of local 
magnetic resonances near the adatom. We show that an 
adatom with inner shell electrons can induce distinctive 
patterns in the nearby LDOS of graphene. This pattern 
can reflect the specific orbital symmetry of the localized 
state, an effect which is unique to the honeycomb lattice. 

In 2 dimensional irreducible representations contain- 
ing doublet states, such as dxy and dx-2-y2 orbitals, we 
find that the LDOS pattern can actually lower the crys- 
tallographic point group symmetry of graphene in the 
situation where the degeneracy of the doublet is lifted by 
interaction effects, which polarize the charge and spin of 
the orbitals^^. In the limiting case where the energy sepa- 
ration of the polarized levels is much larger than the level 
broadening, which is usually the case when the doublet 
state sits near the Dirac points, where the DOS vanishes, 
the orbitals remain orthogonal to each other, and their 
contribution to the LDOS can be detected separately 
by energy resolved STM measurements. In this limit, 
which physically corresponds for instance to the case of 
adatoms in graphene with total spin 1/2, the problem 
can be effectively described by a single orbital Hamil- 
tonian. Although the treatment is generic and can be 
straightforwardly applied to degenerate orbitals in higher 
angular momentum irreducible representations (which in 
graphene include doublet states), to illustrate the main 
effect, we specialize in the case of adatoms with total spin 
~ 1/2, where the adatom induced symmetry patterns in 
the LDOS are more pronounced. 

Furthermore, we discuss the emergence of non-trivial 
particle-hole asymmetries in the energy dependence of 
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the level broadening A(cc;), depending on the particular 
symmetry and position of the localized state in the lat- 
tice. We also address the theory of scanning tunneling 
spectroscopy developed in ref.^^ for the case of 5-wave 
orbit als, and generalize it to describe higher angular mo- 
mentum states in the case of interest, for adatoms with 
total spin near 1/2. 

The outline of the paper is as follows: in section II, 
we describe the generic zero dimensional Hamiltonian of 
an adatom in graphene; in section III we briefly describe 
the role of the orbital symmetry into the formation of 
local magnetic moments and we show the manifestation 
of those orbital symmetries in the LDOS, whenever the 
adatom hybridizes with two or more carbon atoms. In 
section IV we describe the STM tip effects in the LDOS 
and we compute the differential conductance accounting 
for the symmetry of the localized orbitals and their po- 
sition with respect to the sublattices. Finally, in section 
V we present our conclusions. 



II. HAMILTONIAN 

The Hamiltonian of a magnetic adatom in graphene is 
described by a sum of four terms, 

H = Hg^Hf^Hv^Hu, (1) 

where 

H, = -tY,4i^^)K{Il,)^h.c (2) 

is the graphene Hamiltonian in tight-binding, with t ~ 
2.8eV the hoping energy between nearest neighbors sites, 
a (b) is a fermionic annihilation operator in the A {B) 
sublattice, with a =t, i indexing the spin. In momentum 
space. 



(3) 



pa 



where = Yl^=i e*^"^* and ai = x, a2 = — x/2 + >/3^/2, 
and as = -x/2-V^y/2 are the lattice nearest neighbor 
vectors. The second term. 



(4) 



is the Hamiltonian of the localized level with energy eo 
measured from the Dirac point, with m the angular mo- 
mentum projection indexing the different degenerate or- 
bitals contained in a given irreducible representation with 
angular momentum / (for instance, the doublets dxy, 
(i^2_^2, with / = 2 and m = In graphene, due 

to the three fold rotational symmetry of each sublattice, 
the crystaline filed anisotropy lifts the degeneracy of the 
orbitals for different values of |m| < leaving pairs of 
degenerated states (doublets) with angular momentum 
projections ±\m\ with m 7^ 0. 





Figure 1: Honeycomb lattice in graphene, with two distinct 
sublattices (black and white circles) . a) Adatom siting on top 
of a carbon atom on sublattice A, and b) sitting in the center 
of the honeycomb hexagon. Red arrows: nearest neighbor 
vectors. 



The third term in Eq. ([T|) gives the hybridization Hamil- 
tonian. When the adatoms sit on top of a carbon atom as 
in the case of H and F atoms, and also simple molecules^! 
such as NO2, 



Hv = Vyal{0)fm,a^h.c. 



for adsorption on a given site, say on sublattice A. 
Adatoms such as transition metals may instead strongly 
prefer to sit in the hollow site^^ at the center of the 
hexagon in the honeycomb lattice (see Fig. 1). In that 
case, the level is coupled to the graphene bath though 
the hybridization Hamiltonian^^ 



h.c. 



(5) 

where Vx^i {x = a, b) are the hybridization amplitudes 
of the adatom with each of the nearest neighbors carbon 
atoms, which are set by the orbital symmetry of the lo- 
calized state. In momentum space representation, this 
Hamiltonian can be written a a^^i^^ 



where 



t 

a^Y^ ^p 



/m,a + h.c. , (6) 



m,pa 



(m) 



b,p 
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(j) 



(7) 
(8) 



with (j) representing summation over the hybridization 
amplitudes of the adatom with the nearest neighbor car- 
bon atoms on a given sublattice, and TV is the number of 
lattice sites in the extended unit cell of the adatom. The 
discrete sum over momenta can be interchanged by a con- 
tinuous integration, — > A J dk , where A = 2/D'^, 
where D ^ 7eV is the bandwidth. For notational reasons, 
we will set N = 1 from now on. 

For adatoms that sit on top of the carbon atoms, (such 
as hydrogen), 14, p = V and V^^p = 0, for adsorption on 




Figure 2: Illustration of d-wave and /-wave localized orbit als 
for an adatom that sits at the center of the honeycomb lat- 
tice hexagon (H site), a) (1^2 _y2 orbital, that corresponds to 
hybridization amplitudes Vx,i = V and Vx,2 — Vx,3 = —V/2 
with the six nearest neighbor carbon atoms on sublattices 
X — at the vertexes of the hexagon [see Eq. JT]), b) 
dxy orbital corresponding to Vx,i =0 and Vx,2 — —Vx,3 = V. 
c) /x(cc2-3y2) orbital, with hybridization amplitudes Va,i = 
— Vfe,i = V. Adatoms on top carbon sites, and adatoms 
in H or S sites with p-wave, d-wave and out of plane /- 
wave orbitals explicitly break the Csv sublattice symmetry 
in graphene (type / orbitals). s-wave and in-plane /-wave 
orbitals sitting on H or S sites are Csv invariant (type // 
orbitals, see text) 



top of an A site and 14, p = and Vb^p = V for di B 
site. When the adatom sits at the center of the honey- 
comb hexagons (i^-site), the strengths of hybridization 
with the six nearest carbon atoms in the tight-binding 
description depend explicitly on the symmetry of the 
orbital: for example, for 5-wave orbitals, Vx^i = V hy 
symmetry, whereas for in-plane /-wave orbitals, the hy- 
bridization amplitudes are anti-symmetric on the two 
sublattices, Va^i = —Vb^i = V. In the first case (s-wave), 
Va^p = V(j)^ and Vb^p = V(j)p whereas in the second (/- 
wave) Va^p = V0p and Vb^p = —V(j)p. In the case of a 
dx2_y2 orbital Vx^i = V, Vx^2 = Vx^s = whereas 
for a dxy orbital, Vx^i = and Vx^2 = — V^,3 = V, and 
so on, as illustrated in Fig. [2l Other interesting cases 
include for instance substitutional impurities (iS-sites)^^, 
where 14, i = for adatoms sitting on A sites and Vb.i = 
for substitutional adatoms on B sites (see Fig. [3]). A sim- 
ilar description can be for instance applied for adatoms 
sitting on bond sites in between two neighbor carbon 
atoms. 

Finally, the last term in Eq. ((Tj) includes the Coulomb 
energy {U) and exchange energy (J) for the electrons in 
the different orbitals^^, 

HjJ = — ^ ^ ^ ^ {Umm' ~ Jrnrn')fm,af'rn,o-fm',af'rn',(7 
(T m' ^m^m' 

+ ^ Umm' fli^^fm,tfii' ,ifm' (9) 



which can be decomposed at the mean-field level 

into Hu T.m,afm,afmAT^m'Umm'nm'-a + 

J2m'y^miUmm' - Jmm')nm/a], whcrC 



'^m,a — {frn,afm,a) 



(10) 



is the occupation of the orbital with angular momen- 
tum projection m, and spin a. The summation is carried 




Figure 3: Csv invariant orbitals for adatoms sitting in the 
center of the hexagon, on H sites (top), and in substitutional 
(S) sites (bottom). On the left: s-wave orbitals, with zero 
angular momentum (m = 0); on the right: in-plane /-wave 
orbitals (m = ±3). In the two cases, the adatoms hybridize 
equally with the carbon atoms on the same sublattice. 



over the degenerate orbitals in a given irredubile repre- 
sentation, which in graphene correspond to the doublets 
m = ±1^2 1, with \lz\ = 0, 1, . . The mean field interac- 
tion can be absorbed into the definition of the localized 
energy level in Eq. 



H 



f 



(11) 



where 

^m,a = ^0 H~ ^ ^ Urnra''^' 



m', — crH~ ^ ^ {Umm' Jmm'^'^m,' 
m'yi^m 



is the spin dependent renormalized energy of the localized 
states in a given irreducible representation. 



III. LOCAL MAGNETIC MOMENTS 

The formation of local magnetic moments can be ad- 
dressed by the self consistent calculation of the occupa- 
tion for up and down spin states in the different orbitals, 
which follows from integrating the DOS from the bottom 
of the band up to the Fermi level /i,-*^ 

nm,a = -- r dujlmGff^^^uj) (12) 

^ J — oo 



where 



Gff,m,a{^) 



Ejj,^(c^) + zO+] (13) 



is the retarded Green's function of the localized electrons, 
GffAr) = -(r[/(r)/t(0)), and 



(14) 



x,y p 
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is the self-energy of the /-electrons, with x = a^b. G^^ 
are the matrix elements of the retarded Green's func- 
tion of the itinerant electrons in graphene, G^^^p{r) = 

— (T[ap(r)ap(0)]) and so on, which are are defined by 



H 



where 



H 



(15) 



(16) 



is the tight-binding Hamiltonian matrix. More explicitly. 



G°f(p,u;) 



-E- 

q;=± 



where 



I^e(^p)cr^^ -lm(0p)cr^ 



2 

xy 



I0P 



(17) 



(18) 



1 is the identity matrix and and are off-diagonal 



2x2 Pauli matrices, namely <j^^ 



1 and cr?. 



The Green's function of the localized electrons can be 
written more explicitly in the following form: 



G 



R 

ffjm,a 



wZ™^(w) - em,a + iA("')(w) + iO+ 



where 



(19) 



(20) 



gives the quasiparticle residue and A(cj) = —lmT^ff{uj) 
is the level broadening of the localized state^^. 



A(-)(o.) = [y(™)]'^|eL-)|25(o.-a|</.p|), (21) 



which is defined in terms of the generic tight-binding 
phases. 



1 



(m) 



(22) 



where F^^) = max(l/jj^). 

Those phases depend explicitly on the symmetry of the 
localized orbital, which reflect in the relative amplitudes 
of hybridization with the surrounding carbon atoms, and 
also on the relative position of the adatom with respect to 
the sublattices, i.e, if the adatom sits on top of a carbon, 
in the center of the honeycomb hexagon, on a bridge 
site or else in a substitutional site. This formulation is 
completely general and can be easily extended to include 
for instance substitutional impurities in double vacancies. 

In the scenario where the adatom sits on top of a 
carbon atom, the level broadening is given by A(cj) = 



7r]/^p(cj), where p{uj) = \uj\/D'^ is the DOS in graphene 
in the linear portion of the spectrum, and therefore A{uj) 
scales linearly with energy. For adatoms that sit on H 
or S sites, the scaling analysis of the level broadening 
allows a classification in two symmetry groups, depend- 
ing on either if the Csy point group symmetry of the 
honeycomb sublattice is preserved by the adatom or not, 
as previously mentioned in the introduction. As illus- 
trated in Fig. [31 when the electrons hop in and out of an 
adatom sitting on H or S sites, they collect phases which 
give rise to quantum mechanical interference among the 
possible hybridization paths. When the amplitudes of 
hybridization of a localized orbital with the three sur- 
rounding carbons on the same sublattice are identical, in 
which case the Csv point group symmetry of sublattice x 
is preserved, the hopping phases interfere and give rise to 
an anomalous energy scaling of the hybridization, whose 
modulus scales now in the same way as the Kinetic en- 
ergy, \Vx^p\ oc |^p|. In that case, the level broadening 
scales as2^i2^ 



A(cj) ^ 7tV^p{uj) 



(23) 



at low energy, as opposite to the conventional case where 
this interference is frustrated and |V^,p| scales to a con- 
stant near the Dirac points. In the later, A(cj) ex p{uj) 
corresponds to the standard case, whereas in the former 
case the damping is super-linear. The first class of or- 
bit als, which we will refer as type II orbit als, include 
m = and m = 3 angular momentum states, such as in s 
and in-plane /-wave orbitals. The standard "ohmic" class 
(type I) by its turn is described by adatoms on top car- 
bon sites and m = ±1 and m = ±2 angular momentum 
orbitals on H or S sites. To be more concrete, the class of 
type I orbitals is represented by adatoms that sit on top 
of a carbon atom, in which case the orbital symmetry is 
not particularly important, and also by adatoms siting at 
H or S sites with localized orbitals in the Ei{dxz-,dyz) (i-e 
m = ±1) and E2{dxy, dx2-y2) (m = ±2) representations 
of (i-wave orbitals and also f^z^, fyz^, fxyz, fz^x'^-y'^) or- 
bitals in H / S sites. The class of type II orbitals is de- 
scribed by 5, dzz^ fz^ orbitals, where m = 0, and also by 
fxix^-3y^)i and fy(3x^-y2) orbitals (m = ±3) in or 5 
sites. The anomalous scaling of the level broadening in 
Eq. ([231) has been verified explicitly by ab initio methods, 
in particular for the dzz orbital of Co on graphene^^. 

The self-energy of the localized electrons [see Eq. [TH 
can be more explicitly written in the form 

J:ff,M = -uj [Z-Huj) - 1] - iA(-)(a;)|e(D - \^\) . 

(24) 

The density of states of the localized level, Pff^^m,a{^) = 
— l/7rImGjj^^(cj) follows from the substitution of Eq. 
([^ into Eq.' ([T51). 



P//,m,a(^) = - 



A(^)(cj)(9(L>- \UJ\) 



^[ujZ-\uj) -e^,,]2+ [A(-)(, 



(25) 
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Figure 4: Level broadening A(uj) as a function of energy, cj, 
for different orbital symmetries. All energies in units of the 
hopping energy t {V/t = 1/3). a) Type I orbitals for adatoms 
sitting on a top carbon site (light/red solid curve) and for 
dxy-wdive (black line) and d^a^-wave orbitals (dashed line) on 
H sites, b) Type II orbitals on H sites. Black curve: in- 
plane /-wave orbital; light/red curve: s-wave orbital. Inset: 
low energy scaling of the level broadening, A{uj) oc for 
\io\/t < 1 (see text), c) Substititutional s//-wave orbitals 
(black curve) and dxy-wa^ve orbital (light/red curve), on S 
sites. 



Figure 5: Comparison of the LDOS between the two different 
sublattices. Top panels: Energy integrated LDOS around an 
adatom (center) sitting on a top carbon adatom site. Scans for 
a) the opposite sublattice of the impurity and b) for the same 
sublattice. Lower panels: Energy integrated LDOS around 
a localized orbital (center) with 5-wave symmetry, when the 
adatom sits in the center of a honeycomb hexagon (H site), 
c) scans for sublattice A and d) B. The two scans are related 
by a TT-rotation. 



In the linear cone approximation, where the spectrum 
is linearized around the Dirac points, t|^K+p| vp up 
to the cut-off of the band with v ~ 6eVA as the Fermi 
velocity, the level broadening for orbitals of type I is 

A/(cj) = A|cj| (26) 

at low energies, where A = 7r{V/D)'^ is the dimensionless 
hybridization parameter, and 

Z-\uj) = {A/n) In |1 - D^uj^l , (27) 

implying that the quasiparticle residue Z ^ vanishes 
logarithmically at low energy. 

In the case of super-linear damping, for type II orbitals, 
the level broadening scales with the cube of the energy 
within the linear cone approximation, 

Aii{lo) = N,A\Lo\yt\ (28) 

and can be orders of magnitude smaller than in the linear 
case when |eo| <C A/'^ = 1, 2 correspond to the number 
of sublattices the adatom effectively hybridizes {Ns = 2 
for H sites and A^^ = 1 for S sites, as shown in Fig. 
[3]). The quasiparticle residue, Z^r, in this approximation 
is given by 

Z-\uj) = 1 + NsA/{7rt^) [D^ + uj^ In |1 - D^/uj^\] . 

(29) 



In Fig. m we show the energy scaling of the level broad- 
ening for the different orbital symmetries. Particle- hole 
symmetry is preserved for adatoms on top carbon sites, 
where A{uj) follows the DOS, and also for adatoms on 
S sites, which effectively hybridize with only one sub- 
lattice. For adatoms on H sites, which hybridize with 
the two sublattices, particle-hole symmetry is explicitly 
broken in the high energy sector {\uj\ > t) hy the off diag- 
onal matrix elements of the hybridization (at low energy, 
the off diagonal terms average to zero in the momen- 
tum integrals). In particular, (ia,2_^2, dxy-wsive orbitals 
(where Va,i = are strongly damped when the en- 

ergy of the localized state is far above the Dirac point 
(cj > t), but otherwise are weakly damped at negative 
energy states (black solid curve in Fig. [3^). Conversely, 
dzx and (i^^-wave orbitals (where = show the 

opposite trend, in agreement with ab initio calculations 
for Co adatoms in graphene^^^^^. In the same way, s and 
in-plane /-wave orbitals, which couple symmetrically and 
anti-symmetrically with the two sublattices respectively, 
show a strong particle- hole asymmetry at high energies, 
as depicted in Fig. [4)3. In FigHh, we show the level 
broadening for the substitutional case, where particle- 
hole symmetry is restored. In all cases, the peaks at 
\uj\ = t are divergences which are reminiscent of the log- 
arithmic singularity of the DOS around the M point of 
the BZ. 
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Figure 6: LDOS around the adatom (center) at a fixed energy (cj = — 0.2eV) for adatoms on H sites. Tfie top row corresponds 
to a (1^2 _y2 orbital and tfie bottom row to an in-plane /-wave orbital, a) and e): Total LDOS; b) and f): LDOS for the minority 
spin; c) and g) LDOS for the majority spin on sublattice A] d) and h) LDOS for the minority spin on sublattice A. 



1. Spin 1/2 adatoms 

In the situation where the energy of a localized dou- 
blet in a given irreducible representation is close to the 
Dirac points, where the DOS vanishes, the degenerate 
orbitals contained in that representation (say dxy and 
dx'2-y2^ with m = ±2) remain nearly orthogonal to each 
other, since the matrix elements for electronic transitions 
between them, which are mediated by the bath, are very 
small. In that scenario, the lowest energy solution for a 
doublet in the limiting case A ^ involves a maximally 
polarized state where one orbital is fully spin polarized 
and the other is empty (hence with total spin 1/2 and 
total charge of one electron), and the energy separation 
between the orbitals is set by the difference between the 
Coulomb and the exchange energies, U — J, which is typ- 
ically of the order of a few eV^^. More generically, the 
criterion for a maximally polarized solution for a doublet 
state is that the energy separation between the orbitals 
is large compared to the level broadening due to the hy- 
bridization of the orbitals, namely U—J ^ A. When this 
criterion is fulfilled, the electronic transitions between the 
orbitals are suppressed and the effective Hamiltonian of 
the problem maps into the single orbital case described 
in the original work by Anderson^^, where the Coulomb 
Hamiltonian becomes the standard Hubbard U term 

Hu = Uf\f^flh, (30) 

and orbital indexes (m) are suppressed everywhere else 
in Hamiltonian ([T]). The energy of the localized orbital 
becomes 

Hf=Y.^aflU, (31) 



where = eo + Uri-cr gives the energy of the spin polar- 
ized level. 

In the single orbital problem, as in the degenerate case, 
the occupation for up and down spin states can be self- 
consistently calculated at the mean field level from Eq. 
(p!2|) . The emergence of a local magnetic moment follows 
from the appearance of a spin polarized state below the 
Fermi level, say, at energy eo+n^t/, and a virtual (empty) 
state at eo + for the majority spin, with + < 1 
due to the Pauli principle. 

The analysis about the formation of local magnetic 
moments and the zero dimensional phase diagram that 
comes out of the single orbital picture has been discussed 
in detail in ref. 5 for the case of type I orbitals. For type 
II orbitals, the physics is qualitatively similar, except for 
the fact that the formation of a local magnetic moment 
becomes exceedingly easy, even at small [/, due to the 
fact that the broadening of the level can be negligibly 
small when eo/t <C 1. 

For simplicity, in the following, we will analyze the 
single orbital model only, which is valid for orbitals in 
the m = state and also for dublets with total spin 
1/2, when the criterion for the maximal polarization of 
a doublet state, U — J ^ is met. The extension 
to the degenerate multi-orbital case, when the orbitals 
in the doublet state are hybridized, is straightforward. 
The extension may also be applied in the non-degenerate 
multi-orbital case, where the orbitals in the doublet state 
show some degree of charge and spin polarization. In 
standard metals, that translates in the condition U — J > 
ttA^^. From now on, we will drop the orbital indexes m 
and consider only spin polarization effects. 
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A. Local DOS 

The local DOS around the impurity can be computed 
directly from the diagonal matrix elements of the elec- 
tronic Green's function in graphene in the presence of 
the adatom, 

^,(r,a;) = -llm ^ ^ G« ,,(p, p', a;)e*(P-P')-^ , 

cr p,p' 

(32) 

where 

GxyAV. p', = ^v.v'G^y^V) + Ax{p)Gff^a{ioj)Ay{p') , 

(33) 

and 

A,(P) = J2 Glyip)Vy,p (34) 

y=a,b 

hip) = J2 Kp^Up)^ (35) 

y=a,b 

with Va^p and Vb^p defined in Eq. ([7j) and ([H]). 

In Fig. [5] we show the topography maps around the 
impurity, which describe the local DOS integrated in en- 
ergy. We use the set of parameters V = leV, U = leV, 
jii = O.leV and eo = — 0.5eV, which correspond to a spin 
polarized state with ^ 0.04 and ^ 0.96. In the case 
where the adatoms sit on top of a carbon site on a given 
sublattice (top panels), the maps show a clear asymmetry 
between the integrated LDOS of the two different sub- 
lattices. The pattern on the opposite sublattice of the 
impurity (Fig. [Sti) has a lower point group symmetry 
than in the same sublattice (Fig. [SJd), what comes from 
the fact that the adatom in this case has only three near- 
est neighbor carbon sites but six next-nearest neighbor 
ones. For adatoms sitting in the center of the honeycomb 
hexagon (lower panels), there is no distinction between 
the patterns of the two different sublattices, except for 
a rotation of tt. Fig. and d depict the integrated 
LDOS for an 5-wave orbital (m = 0) sitting on an H 
site. The intensity of the integrated LDOS maps is also 
much weaker in the lower panels compared to the upper 
ones, reflecting the fact that the hybridization for H or S 
sites is mediated by hopping, and hence weaker than in 
the top carbon site case for the same set of parameters. 

The analysis of the LDOS also permits to identify the 
symmetry of the localized orbital when the adatom sits 
either on H or S sites. In Fig. [Hlwe show the LDOS at 
fixed energy for both a d^'^^y'i orbital (Fig. [6^-d) and a 
fx{x'^-3y'^) state (Fig. [Ht-h). In the former, the orbital 
(i-wave symmetry of the localized state has a clear fin- 
gerprint in the induced DOS nearby the adatom. The 
signature is specially pronounced when the LDOS is re- 
solved for the minority spins, as noticed by direct com- 
parison of Fig. [6)3 and f. In panels c), d) and g), h), 
we show the distinct patterns for the LDOS of majority 
and minority spins on a given sublattice. In the case of 
an in-plane /-wave state, which explicitly preserves the 



Top -sublattice B Top -sublattice A 




H site (d-wave) H site (f-wave) 




Figure 7: Fourier transform of the energy integrated LDOS 
around the adatom. Solid hexagon line indicates the Brillouin 
zone. The peaks at the center of the zone correspond to for- 
ward scattering processes, whereas the peaks at the corners 
of the zone (K) points correspond to backscattering between 
the two valleys. Top panels: adatom on a top carbon site 
case; Fourier transform of the LDOS for the a) opposite and 
b) same sublattice of the adatom. c) superposition of the 
patterns in panels a) and b), for both sublattices. d) s-wave 
orbital at an H site. Lower panels: e) (ia,2_^2-wave orbital 
and f) /cc(cc2_3^2) orbital also at H sites. The insets in panels 
d) , e) and f) show the details of the forward scattering peaks 
at the center of the BZ. 

point group symmetry of the sublattices, the fingerprint, 
although more subtle, can be easily characterized by a 
Fourier analysis of the LDOS, which maps the scattering 
wavevectors responsible for the emergence of local reso- 
nances nearby the adatom. 

In Fig. [3 we plot the corresponding maps of the 
energy integrated LDOS in the reciprocal space. The 
solid hexagonal line indicates the Brillouin zone (BZ). 
In Fig. [7^^, we show the maps for the opposite sub- 
lattice of the impurity, whereas on Fig. [JJd, we depict 
the Fourier transformed maps for the same sublattice 
of the adatom. The central peaks around the center of 
the zone (F point) indicate forward scattering processes 
that contribute to the resonant states nearby the adatom, 
whereas the peaks centered at the K points at the edges 
the BZ indicate backscattering processes, which connect 
the different valleys. For the opposite sublattice of the 
adatom (Fig. [7^), the backscattering processes at the 



K point and the forward scattering ones at F are sig- 
nificantly attenuated. In the reciprocal space maps for 
the same sublattice of the impurity (Fig. [Jb), where 
unitary scattering should dominate, backscattering pro- 
cesses are strongly enhanced, followed by the presence of 
subdominant forward scattering peaks. In panel [7h, we 
depict the Fourier transformed map for both sublattices. 
In those plots (top carbon case), the amount of scatter- 
ing at the M points, which indicate the position of the 
Van Hove singularities, is weak compared to the other 
dominant processes. The shape of the forward scattering 
peak at the F point also reflects the symmetry of the hy- 
bridization matrix elements in the Hamiltonian. In the 
top carbon case, the F peak is isotropic. 

In Fig. [TJi we show the reciprocal space maps for the 
energy integrated LDOS for an 5-wave orbital siting at 
an H site. In this case, the height of the central peak 
is significantly small compared to the dominant peaks 
around the K points, indicating strong enhancement of 
the backscattering compared to forward scattering pro- 
cesses. In panel [7^ we depict the case of a dx2_y2-wSiYe 
orbital also at an H site, whereas in panel [IF we show the 
signature of an fx{x^-3y^) orbital {H site) in the recipro- 
cal space. For d and /-wave orbit als, destructive interfer- 
ence leads to attenuation of the backscattering peaks at 
the K points, in particular in the d-wave case. In the in- 
sets of Fig. [7]i and f we show in detail the features of the 
forward scattering peaks for 5-wave and /a?(a;2_3^2)-wave 
orbitals, respectively. Both peaks reflect the underlying 
Csy symmetry of the sublattices, which are incorporated 
into the hybridization matrix elements of the Hamilto- 
nian for orbitals of type II. For a (i-wave orbital (inset of 
Fig. C^), the forward scattering peak has C2v symmetry. 



IV. STM TIP EFFECTS 

Let us now consider a problem of more practical inter- 
est for tunneling microscopy experiments, where we addi- 
tionally include an STM tip close to a spin 1/2 impurity. 
The electrons in the metallic tip follow the Hamiltonian 



(36) 



where ep = /2m is the electronic dispersion of an elec- 
tron gas, with m the effective mass. The electrons can 
tunnel either to the carbon sites in graphene or to the 
impurity. In the former case, the tunneling process is 
described by the Hamiltonian 

Hg-t = Y,Y.taal{Ki)c,{Ki - r) 

(i) a 




Figure 8: Schematic drawing of the STM tip nearby an 
adatom (small light [red] circle) on top of graphene. Black 
and white circles: carbon atoms on sublattices A and B. is 
the in-plane distance of the impurity to the tip and z the out 
of plane distance from the center of the tip to the graphene 
layer. 



over the tip nearest neighbor carbon sites (R^) on a given 
sublattice, and r = (i?, z) is the position of the center of 
the tip, where R is the horizontal distance of the tip 
to the impurity and z is the distance of the tip to the 
graphene layer. 

The single-particle wave- functions t/jt^pir) describing 
the electronic state at the tip, namely 

c^(Ri - r) = ^ V^t,p(Ri - r)c^p , 
p 

can be expanded in spherical waves from the center of the 
tip, '0^,p(r) (X e~'^p^/r, where c^p is a second quantized 
operator for the tip electrons. The factor— 



2m{(j)t 



3) 



(38) 



gives the effective tunneling barrier between the tip and 
the rest of the system, and is defined by the electronic 
work function of the tip, (j)f Since the single-particle 
wave functions of the graphene electrons can be expanded 
in plane waves, ?/;^,k(Ri) = e*^ -^*^^, the Hamiltonian 
(|57|) becomes 



Hg-t = XI [^«'kp(r) al^c^p + t6,kp(r) bl^c^p 

(J kp 



where 



t.,kp(r) ~ "JHe-^-e-'''-^ = 



^X,p(^) 6 



-ik R 



+ /l.C. , 

(39) 
(40) 



describes the spatially averaged hopping matrix elements 
between the tip and graphene, where the position of each 
of the carbon atoms underneath the tip is effectively re- 
placed by the in-plane position of the center of the tip 
with respect to the impurity. 

The Hamiltonian for the tunneling from the tip to the 
impurity is given by: 



where ta^h are the electronic tunneling energy from the 
tip to sublattices A, and B in graphene, {i)t denotes sum 



Hf.t=tfY,flc.{r) + h.c. 

a 



(41) 
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where t/_t is the tunnehng energy from the tip to the 
impurity. In a similar way, we can write 

Hf_t = ^tj,p{r)fXp + h.c., (42) 

where 

*/,pW=^/^e-«pl-l. (43) 



where |r| = \^ R'^ -\- z'^ measures the center of the tip with 
respect to the position of the impurity. We will assume 
that the tip is sufficiently large such that local gating 
effects due to the proximity of the tip to graphene can 
be effectively absorbed into the local definition of the 
chemical potential nearby the adatom^^. 

As a brief comment about orders of magnitude for the 
several quantities, in most materials, the typical work 
function (j)t is of the order of a few eV. In the case where 
the effective mass m is of the same order of the bare mass 
of the electron, z^"^, translates into a characteristic tun- 
neling length scale topically larger than 1 nm. Also, since 
STM tips have a typical radius of the order of Inm, 
which accounts for the distance between the center of the 
tip to graphene is typically a number of the same order. 
A more detailed description about the characterization 
of the tip and the surface in the STM problem can be 
found in Ref. il9i . 

A. Green's functions 



the presence of the STM tip. A detailed calculation of 
the equilibrium and also out of equilibrium Green's func- 
tions in the STM problem for metallic surfaces can be 
found in Ref. 

To further simplify matters, we assume here that A^p ~ 
Hi is computed at the Fermi energy and hence is momen- 
tum independent, in which case 

t:c,kp(r) ^ U{z) e-^^-^ = t,,k(r) , (45) 

with X = a^h. Since the STM tip is typically large com- 
pared to the lattice spacing in graphene, one can further 
simplify things by assuming p = ^6 p. We will keep the 

0.54 0.58 0.62 '^.5 0.6 




Sublattice B Sublattice A 




0.06 0.11 0.16 0.06 0.11 0.16 



Now we generalize the original Hamiltonian of the An- 
derson problem to include Hamiltonian terms (|36|). ([39]) 
and (|12D, 

H = Hg^Hf^Hv^Hu^Ht^ Hg^t + Hf^t • (44) 

In the following, we will assume perturbation theory in 
the hybridization of the tip with the rest of the system, 
namely tf and tx {x = a, b) are small compared to the 
hybridization V of the adatom with the host (graphene). 
This is not nevertheless a strict requirement at the mean 
field level, and the 'exact' expressions of the Green's func- 
tions are shown in the Appendix. In addition, we will 
also assume that the system remains in equilibrium in 



Figure 9: Energy integrated LDOS around the adatom (cen- 
ter) in the presence of an STM tip. a) Scans for the same sub- 
lattice of the impurity and b) for the opposite sublattice (top 
carbon site case). Scans for sublattice c) A and d) B, nearby 
an adatom sitting at the center of a honeycomb hexagon for 
some magnetic m = angular momentum state (center). 



a, b labels below for completeness. A similar assumption 
will be made for the tunneling matrix element between 
the tip and the adatom, tf^p{r) tfi^)- 

The matrix elements of the renormalized Green's func- 
tion is the a, b sublattice basis can be calculated straight- 
forwardly. 



GxyAP. P', = ^P,P'Gly{p) + r,(p, r)Ax{p)GffAiuj)Ay{p')fy{p',r) + T,(p, r)fy{p\ r) ^ G^ik) , (46) 

k 



The quantities A^; (p) and A^^ (p) were defined in Eq. and (|35]) , whereas 

T,(p,r) = G%ip)ty,^ir) (47) 

y=a,b 

f,ip,r) = t;jr)Glip) (48) 

y^a,b 
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contain the tunneling amplitudes and phases for the elec- 
trons as they hop between the tip and the A, B sublat- 
tices. G^^(k,r) = -(T[ck(r)4(0)]) is the bare Green's 
function of the electrons of the tip, 



lUJ — Cr- 



(49) 



while r and its conjugate form F in Eq. define the 
vertex renormalization due to the presence of the tip, 

r.(p,r) = l + ^^t^(r,a;)^GL(fc) (50) 

f,(p,r) = l + TdP2^tf{r,u;)Y,Gl{k), (51) 

where the quantities 

tf{r,iu) = tf{r)+ J2 Ay{-'R,i^)ty{z) (52) 



y^a,b 



if{r,iuj) = tf{r) + t*y{z)Ky{Ii,iuj), 

y=a,b 



(53) 



give the renormalized tunneling functions between the 
tip and the adatom, whose bare form, tf{r)^ is defined in 
Eq. (gSl) . We also defined 

A,(R) = ^A(fe)e^^-^ 

k 

A,(R) = ^A(fc)e*-^ 



which refiects the infiuence of the tip into the wavefunc- 
tion of the localized states. 

The imaginary part of the self energy, ImS^j (r, zcj), 
renormalizes the level broadening A(cj), defined in Eq. 
dnH), due to the hybridization of the localized electrons 
with the electrons in the tip. In contrast with metallic 
hosts, which have a large DOS, in graphene the metal- 
lic tip can locally overwhelm the hybridization of the 
adatom with the nearby carbon atoms. In the situation 
where the level broadening becomes large enough as to 
overcome local correlation effects in the localized state, 
the tip might eventually lead to suppression of the local 
magnetism. This effect will be discussed in more detail 
in sec. IV. B. 

Finally, other useful quantities are two off diagonal 
Green's functions Gcx^aiPi'^) = ~(^[c(t)x'''(0)]), with 
X = a^b, which are given by 



(57) 



and also 



GcfAp) = G°e(p)i/(r, iu)GffAi^) , (58) 

which are required for computing the differential conduc- 
tance (see sec. IV. B.l). 



as the Fourier transforms of A(/c) and A(/c) [see Eq. 
and 

The self-energy correction to the localized electrons, 
i;//(cj), as given in Eq. ([Tlj), is also dressed by the prox- 
imity of the STM tip and assumes the form 



^// (r, = ^ff {i^) + S// (r, iuj) 



(54) 



where 



(r, ico) = tf{r)ifir) ^ G,,{k) . (55) 



gives the contribution from the tip to leading order in 
tf and tx. In the presence of the STM tip, the Green's 
function of the localized electrons depends explicitly on 
the distance between the tip to the adatom, 

= [^^ -^cr- S//(r, iuj)]~^ , (56) 



L Local DOS 



Besides the localized state, the local DOS around the 
adatom is also affected by the presence of the STM tip. 
The local DOS nearby the impurity is also indirectly af- 
fected by the hybridization of the orbitals of the tip with 
the adatom localized orbital. For instance, for a mag- 
netic adatom, the DOS is expected to be spin polarized 
on a given sublattice x, 

p^Ar, ^) = --Im ^ G^^p, p', u;)e^(P-P')-^ , (59) 

P,P' 



where the diagonal Green's function Gxx,a{p^p' ^'f^^) is 
explicitly shown in Eq. In a more explicit form. 



p^,^(r,cj) --Im 

TT 



J2 GlM + [r. * AJ(r) GffA^) [f , * AJ(r) + T,{z)%(z) ^ Gl,(k) 



(60) 
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Figure 10: Graphene LDOS at the adatom site (top carbon 
case) for tf/ta — 0.1, 0.7, 1.1 and 1.6, from (a) to (d) (ta = 
th = 0.15eV). Black curve: n^; brown: n^. Total LDOS: 
+ n^. The two resonant peaks have energy eo -\- n\U and 
eo + n^U ^ with eo = —0.5 eV and U = leV (see text). 



gives the local DOS per spin, where 

[r,*AJ(r) = ^e'P-^r,(p,r)A,(p) (61) 

P 

[f,*AJ(r) = ^e-*P-^f,(p,r)A,(p) (62) 



state. For moderate [/, this effect may imply in the even- 
tual suppression of the localized moment when the tip is 
sufficiently close to the adatom, permitting the detection 
of the local moment with a non-magnetic tip away from 
the Kondo regime. This effect has also a clear manifesta- 
tion in the differential conductance curves of STM, which 
we describe in the next subsection. 

We remark that in the regime where the STM tip leads 
to a suppresion of the local moment, the single orbital 
picture we considered so far breaks down for orbitals 
with \m\ > 1 when the energy splitting of the spin po- 
larized dublet states in a given representation (say dxy 
and dx2_y2) becomes comparable to the level broadening 
induced by the tip. In that regime, the suppression of 
the local moment is followed by a strong enhancement of 
the hybridization among the orbitals, and a multi-orbital 
description is needed. This extension with nearly degen- 
erate orbitals is straightforward and will be considered 
somewhere else. 



B. Tunneling current 



is the Fourier transform convoluted over the the product 
of the r(j9, r) and Ax{p) functions (and their respective 
conjugate forms), as defined in Eq. (|35|) and ([5Q|) . 

(|5T]) , while Tx{z,uj) is by definition T(r, cc;)|j^^q [see Eq. 
dlTj)], and hence independent of the horizontal distance 
between the tip and the adatom. 

In Fig. we show the topography maps for LDOS 
in the presence of the STM tip for both sublattices in 
the case of an adatom sitting on top of a carbon atom 
(top panels) and also for an adatom on an H site (lower 
panels). In those plots, we use the same set of parameters 
as before, V = leV, U = leV, /i = O.leV, eo = -0.5eV, 
and additionally the parameters = 4 eV for the band 
width of the tip and cd = 2eV for the Fermi energy 
of the tip. The tunneling parameters between the tip 
and the system where chosen to be = 0.02 eV and 
ta = tb = 0.2eV. As expected, the plots show basically 
the same qualitative features as the ones shown in Fig. 
[5] for the actual DOS on graphene in the absence of the 
STM tip. 

The most important feature of the tip in the LDOS 
can be observed in the evolution of the magnetic reso- 
nances as the STM tip gets close to the adatom. In Fig. 
nm we show the LDOS measured at the position of an 
adatom (R = 0) that sits on top of a carbon atom. Us- 
ing the same set of parameters as before, we show the 
evolution of the two spin polarized resonant peaks at en- 
ergy eo-\-n^U and eo + n|/7 as we increase the strength of 
the hopping amplitude between the tip and the adatom, 
tf. As tf increases (panels [TDb to d), the peaks are pro- 
gressively suppressed as the high energy peak for virtual 
n^states is red shifted. As mentioned in the previous sec- 
tion, this effect stems from the fact that due to the low 
DOS in graphene, the proximity of the metallic STM tip 
can substantially enhance the broadening of the localized 



The tunneling current from the tip is defined by^^: 



I = -e 



dt 



(63) 



where Nc = J^ka ^kcr^ka is the number operator for the 
c electrons in the tip, and e is the electron charge. The 
motion equation for this operator is 



dtNc 



where H is the full Hamiltonian defined in Eq. in- 
cluding hopping matrix elements between the tip and the 
system. After a straightforward algebra, the total current 
follows from the sum of three different contributions that 
arise from tip tunneling processes to either sublattices A, 
and B or else to the adatom localized state, 



/ = -2elm 



f 

J — C 



2^ 



k,cr 



kp,t7 x—a,b 



where 



G>{t,t') 



,(64) 



(65) 
(66) 



are real time "lesser" and "greater" Green's functions, 
which should be distinguished from retarded (G^) and 
advanced (G^) ones, which are time ordered. If A{ij) = 
B{uj)C{uj)^ one may show that^^ 

A<{uj) = B^{uj)C<{uj) + B<{uj)C^{uj) , (67) 
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in which case the total current can be written as 



/ = — 2elm / X 

27rz 



'f 

J — C 



J2 [G'J{k)J:<{r,uj)+Gl<{k)^i{r,u)] (68) 



where 



p x=a,b 



+tf{r)tf{v)J2GffAi^) (69) 



is the self energy correction to the Green's function of 
the c-electrons in the tip, 

Gcc(p, p', i^) = [{i^ - ep)(^p,p' - Scc(r, iuj)]~^ . (70) 

The first term in Eq. (|69|) gives the self-energy contribu- 
tion due to the graphene electrons, while the second one 
is the contribution from the adatom. Using the fluctua- 
tion dissipation theorem^^ 



G<{co)=if{oj)A{co)., 



(71) 



where A{ijj) = —21mG^{uj) is the spectral function, the 
total current is given by 

/CO 
dupt{r,uj)pe{uj') [f{u)-f{oj')], 
-co 

(72) 

where f{uj) = [e^/^ + 1] ^ is the Fermi distribution, pc 
is the DOS at the STM tip. 



and pt is defined as 

P*(r,w) = 



ImSf,(r,o.). 



(73) 



(74) 



This term has units of DOS, and accounts for the phases 
acquired by the electrons in the tunneling process be- 
tween the tip and the localized state of the adatom. It 
can be conveniently rewritten in the following form: 

{TT 

^ [(77 - m) ImG'//,a + (^7 + ^7)ReG//,a] | , 



where 



7rFpo(^ 



^Im^ T,(R = 0,c^). (75) 



Po(^) is the bare local DOS of graphene in the absence 
of the impurity and the tip, and V = max(14^i), with 



Vx^i the hybridization amplitudes of the adatom with the 
nearest carbon atoms, as defined in Eq. ([22) • The other 
parameter, q and its conjugate form, q are the Fano fac- 
tors. 



q{r,uj) = 



1 Ret/(r, c 



ta{z)V 7rpo{uj) 



(76) 



while 



-/ N 1 Retf(r, cj) 

ta{z)V 7Tpo{oJ) 



7(r,c^) = - . ImA,(R,c^) (78) 

7(r,c^) = - . Yl ImA,(R,c^) (79) 



x=a,b 



gives the corresponding damping factor. These factors 
characterize the Fano resonances in the differential con- 
ductance curves in the vicinity of a localized state. 



1. Differential conductance 

The Green's function of the localized electrons can be 
written in a more compact form as: 



(80) 



where S//(r, uj) is the dressed self-energy of the localized 
electrons due to the proximity of the tip, as defined in 
Eq. ([54j) . and ^cr(^) is defined as 



^a(r,cj) 



Cg - Rei://(r,cj) 
ImS//(r,a;) 



(81) 



The differential conductance follows by computing I /Vb 
in the limit of Vb = oj' — OJ ^ 0. Since d/(co'')/dco' = 
—6{uj') at zero temperature, the differential conductance 
can be written in the more standard form 



27repc{0)t}{z)po{uJb) x (82) 
"gg- 77+ ((77 + ^7)^^ 



where ujb is the bias voltage. The first term in parenthesis 
defines the DC due to the DOS in graphene. The second 
one is explicitly defined in terms of the Fano parameters 
and gives the contribution due to the presence of the 
magnetic adatom. 

The experimental detection of a localized state with 
STM tips is based on the principle of quantum interfer- 
ence between the two different hybridization paths the 
electrons can take when they tunnel from the impurity 
to the localized state. In one way, the electrons can tun- 
nel directly to the localized state. On the other, they 



13 



can also tunnel to the host material (graphene) and then 
hybridize with the localized orbital. The signature of 
such interference appears in the differential conductance 
curves in the form of a Fano resonance. In graphene, 
the electrons have additional sublattice quantum num- 
bers which may give rise to additional interference ef- 
fects, depending on the position of the adatom relative 
to the two different sublattices. In the case where the 
adatom sits in the center of the honeycomb hexagon, for 
a given sublattice, there are three different paths the elec- 
trons in graphene can take to hybridize with the adatom. 
Destructive interference between the different paths in a 
given sublattice can suppress the Fano character of the 
resonance and change the shape of the DC curves. 

When the tip is above the adatom (i^ = 0), the con- 
jugate forms q = q and 7 = 7 in Eq. ([76 |) - ([79|) are the 
same. In the simplest scenario, where an adatom sits on 
top of a carbon atom, say on site A, the Fano factor is 
defined explicitly in terms of the self-energy for orbitals 
of type I, Rei;jj(cj) = uj[ZY'^{oj) — 1], as given in Eq. 
([57j), namely 



Type I orbitals 



Type II orbitals 



Ve + {ta{z)/V)ReEjf{i0) 

nta{z)Vpo{uj) 



(83) 



The damping in this case is 7^ = 1, by noticing that in- 
tegrals with off diagonal matrix elements of the Green's 
function, such as Gab = 0. In a different scenario, 
for adatoms of type II, which sit either in 5* or sites 
and possess orbitals with Csy point group symmetry, 
as discussed in Sec. Ill, the hybridization matrix ele- 



ments have the form V^^p = ±F^p and Vc 



a,p 



in which case one can easily check that A(R = 0, w) 
'Ex,y=a,b Ek GxyVy^p = 0. In that case, 



(84) 



and 7^^ = 0. In the more generic case, for type I orbitals 
(the ones which are not C^v invariant) that sit on H 
or S sites, the damping factor 7 interpolates between 
and Ns = 1,2, the number of sublattices the adatom 
effectively hybridizes. 

As the usual theory of Fano resonances^^, when q/^ ^ 
1, the DC curve has the form of a peak, whereas in the 
opposite limit, for ^^77 <C 1 it is shows a dip. In Fig. [TT] 
we show the DC induced by the presence of the adatom 
in the case where it sits on top of a carbon site (left 
panels) and also for an type II orbital at an H site (right 
panels). In panels a) and b) we assume a fixed set of 
parameters as in FigfTOl and change the ratios tf/ta for 
a given fixed value of ta = 0.1 5eV. For tf/ta ^ 0.1, the 
DC curve of type I orbitals has a Fano shape, whereas 
for type II orbitals the Fano resonance is suppressed (see 
insets of Fig. [TT]) . 

The evolution of the separation of the peaks follows the 
same trend shown before in the LDOS (see Fig. [T0|). The 
increase of tf leads to a gradual suppression of the local 
magnetic moment, and as a consequence to a decrease in 
the separation of the two peaks. When is large enough. 
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Figure 11: Differential conductance induced by the adatom 
versus bias when the tip is right above the adatom (R = 0). 
The different curves are shifted vertically. Left: type one 
orbitals, for adatoms that sit on top of carbon site, or else 
in if or 5* sites, for orbitals that break the Csv symmetry of 
the crystal; (right) Csv invariant orbitals at 5* or if sites. See 
details in the text, a), b) ta — th — 0.15eV and tf/ta — 1-6, 
1.5, 1.35, 1.1, 0.7, 0.25 and 0.1 (inset), from top to bottom. 



the two DC magnetic peaks eventually merge on top of 
each other, destroying the local moment completely. In 
this situation, the level broadening of the localized state 
is controlled by the proximity of the tip, and the local 
magnetic moment can be in principle turned off by bring- 
ing the tip closer to the adatom. Additional fluctuations, 
which are expected to play a role at very large t / , should 
further enhance the suppression of the local moment by 
the STM tip. This effect is inhibited by strong corre- 
lations, at large [/, which make the local moment more 
resilient to the influence of the tip. Gating effects that 
bring the Fermi level close to the energy of the local- 
ized state may on the other hand favor a mixed valence 
regime state, where the local moment can be more easily 
destroyed by the tip even at very large U . 



V. CONCLUSION 

Unlike the case of metallic hosts, in graphene the sym- 
metry of the localized orbital has clear fingerprints in 
the LDOS nearby the adatom whenever the adatom, hy- 
bridizes with two or more carbon atoms. We showed that 
the real and momentum space STM scanning maps can 
reveal not only the position of the adatom with respect 
to the sublattices but can also indicate the orbital sym- 
metry of the localized state and possibly its magnetic 
state. 

We have described in detail how sublattice quantum 
numbers in combination with orbital symmetry effects 
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influence the Fano resonances in the differential conduc- 
tance nearby the adatom. To ihustrate the effect, we 
analyzed the problem in the single orbital picture, which 
is valid for orbitals in one dimensional irreducible rep- 
resentations (such as s-wave, d^^wave, etc), as well as 
in the more interesting case of 2 dimensional irreducible 
representations, such as in the doublets {dx2_y2^ dxy), 
(/x(a;2-3y2),4(3y2_^2)), etc, when the energy separation 
of the orbitals in the doublet state, driven for instance 
by spin polarization effects, is large compared to the level 
broadening. The formalism is generic for any orbital sym- 
metry state, and can be trivially extended to describe the 
degenerate case as well. 
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Appendix A: Exact Green's functions 

At the mean fleld level, the Green's functions described 
in sec. IV can be written in an exact form by solving the 
equations of motion for the fermionic operators ak,cr, ^k,cr, 
Cp^a and fa. After a cumbersome but straightforward 
algebra, the flnal expressions are. 



G:c2y,a(p,P^^cJ) =^p,p.G^^(p)+r^(p,r)A^(p)G//,,M (Al) 

(A2) 



Gcx,a{p^p\i^) = ^Gcc(p,k',zcj) [T^,k(r,p ) + t/,kr:c,p'G//,^(icj)A^(p ) 
k 

I 



and 



where T>ff{r,i(jj) = T>ff{iuj) + Iljj{r, iuj) is the corre- 
sponding self-energy, with 



(A3) ^ff{r,iu) = ^t/,k(r)Gcc(k,k^^c^)t7k'(r) (A9) 



kk' 



where p = (p^iuu). The quantities Ax{p) and Ax{p) have 

the same deflnitions as before [see Eq. (|3ll) and (|35], contribution of the tip. The other quantities in- 



whereas 



Tp,,,(r,p) = ^ G°^(p)t,,p,p,(r) 



(A4) 



y=a,b 



rp.,,(r,p) = ^ t;p,p,(r)G°,(p). (A5) 

y=a,b 

The Green's function of the c-electrons is deflned as 

Gcc(p, p', i^) = [{i^ - ep)^^pp' - Scc,pp'(^, ^^)]""^ 1 

(A6) 

where Sec is the self-energy for the electrons in the tip due 
to hybridization effects with the electrons in graphene 
only. 



elude the renormalized hybridization of the tip with the 
adatom, 



k y=a,b 

i7,p(r, iu:) = tf,p{r) + ^ ^ t^kp^k) , (All) 

k y—a^b 

and 



^xip) = 1 + ^t/,k(r,zcj)Gcc(k,k',za;)^ 



kk' 



^ k^x(r,p) 



k x—ttjb 

The Green's function of the /-electrons is given by 

GffA^^) = [^^ -^a- ^ff{r, ic^)r^ , (A8) 



1 ^x{p) 



kk 



which are vertex corrections that appear in Eq. (|Aip and 
T2I). 
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